Effects of four-body breakup on 6 Li elastic scattering near the Coulomb barrier 
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We investigate projectile breakup effects on 6 Li+ 209 Bi elastic scattering near the Coulomb barrier with 
the four-body version of the continuum-discretized coupled-channels method (four-body CDCC). This is the 
first application of four-body CDCC to 6 Li elastic scattering. The elastic scattering is well described by the 
p+n+ 4 He+ 209 Bi four-body model. We propose a reasonable three-body model for describing the four-body 
scattering, clarifying four-body dynamics of the elastic scattering. 
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Introduction. Plenty of nuclei are considered to have two- 
cluster or three-cluster configurations as their main compo- 
nents. Three-cluster dynamics is, however, nontrivial com- 
pared with two-cluster dynamics. Systematic understanding 
of three-cluster dynamics is hence important. There are many 
nuclei that can be described by three-cluster models. For ex- 
ample, low-lying states of 6 He and 6 Li are explained by N 
+ N + 4 He three-body models lUHfjl], where N stands for a 
nucleon. The comparison of the two nuclei is important to 
see the difference between dineutron and proton-neutron cor- 
relations. Two-neutron halo nuclei such as n Li, 1 Be, and 
22 C are reasonably described by an n + n + X three-cluster 
model, where X is a core nucleus. Properties of these three- 
cluster configurations should be confirmed by measuring scat- 
tering of the nuclei and analyzing the measured cross sections 
with accurate reaction theories. The reactions are essentially 
four-body scattering composed of three constituents of pro- 
jectile and a target nucleus. Accurate theoretical description 
of four-body scattering is thus an important subject in nuclear 
physics. 

The continuum-discretized coupled-channels method 
(CDCC) is a fully quantum-mechanical method of describing 
not only three-body scattering but also four-body scatter- 
ing 0-121 • CDCC has succeeded in reproducing experimental 
data on both three- and four-body scattering. The theoretical 
foundation of CDCC is shown with the distorted Faddeev 
equation llTol - [l2tl . CDCC for four-body (three-body) scat- 
tering is often called four-body (three-body) CDCC; see 
Refs |[l3l - l25ll and references therein for four-body CDCC. So 
far four-body CDCC was applied to only 6 He scattering. 

For 6 He + 209 Bi scattering at 19 and 22.5 MeV near the 
Coulomb barrier, the measured total reaction cross sections 
are largely enhanced in comparison with that for 6 Li + 209 Bi 
scatterin g a t 29.9 and 32.8 MeV near the Coulomb bar- 
rier lUSlSll. Keeley et al. JH analyzed the 6 He + 209 Bi 
scattering with three-body CDCC in which the 6 He + 209 Bi 
system was assumed to be a 2 n + 4 He + 209 Bi three-body 
system, i.e., a pair of extra neutrons in 6 He was treated as a 



single particle, dineutron ( 2 n). The enhancement of the total 
reaction cross section of the 6 He + 209 Bi scattering is found 
to be due to the electric dipole (£1) excitation of 6 He to its 
continuum states l29ll . i.e., Coulomb breakup of 6 He, which 
is almost absent in the 6 Li + 209 Bi scattering. The three-body 
CDCC calculation, however, does not reproduce the angular 
distribution of the measured elastic cross section and overesti- 
mates the measured total reaction cross section by a factor of 
2.5. This problem is solved by four-body CDCC 1U9I1 in which 
the total system is assumed to be an n + n + 4 He + 209 Bi four- 
body system. 

The 6 Li + 209 Bi scattering near the Coulomb barrier was, 
meanwhile, analyzed with three-body CDCC by assuming a 
d + 4 He + 209 Bi three-body model JH]. The three-body 
CDCC calculation could not reproduce the data without nor- 
malization factors for the potentials between 6 Li and 209 Bi. 
This result indicates that four-body CDCC should be applied 
to the 6 Li + 209 Bi scattering. 

In this paper, we analyze 6 Li + 209 Bi elastic scattering at 
29.9 and 32.8 MeV with four-body CDCC by assuming the 
p + n + 4 He + 209 Bi four-body model. This is the first ap- 
plication of four-body CDCC to 6 Li scattering. The four- 
body CDCC calculation reproduces the measured elastic cross 
sections, whereas the previous three-body CDCC calculation 
does not. Four-body dynamics of the elastic scattering is in- 
vestigated, and it is discussed what causes the failure of the 
previous three-body CDCC calculation. Finally we propose 
a reasonable three-body model for describing the four-body 
scattering. 

Theoretical framework. One of the most natural frame- 
works to describe 6 Li + 209 Bi scattering is the p + n + 4 He 
+ 209 Bi four-body model. Dynamics of the scattering is gov- 
erned by the Schrodinger equation 

{H-E)V = (1) 

for the total wave function '!?, where E is a total energy of the 
system. The total Hamiltonian H is defined by 

H = K R + U + h (2) 
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with 



U = U n (R n ) + U P (R P ) + U a (R a ) + 



e 2 ZuZBi 
R : 



(3) 



2 



where h denotes the internal Hamiltonian of 6 Li, R is the 
center-of-mass coordinate of 6 Li relative to 209 Bi, Kr stands 
for the kinetic energy operator associated with R, and U x de- 
scribes the nuclear part of the optical potential between x and 
209 Bi as a function of the relative coordinate R x . As U a , we 
adopt the optical potential of Barnett and Lilley ll30ll . Param- 
eters of U n are fitted to reproduce experimental data ll3lTl on 
n + 209 Bi elastic scattering at 5 MeV, where only the central 
interaction is taken for simplicity. As shown in Fig. Q] the 
neutron optical potential U® p thus fitted is consistent with 
the data. The resultant parameter set is the same as that in 
the global optical potential of Koning and Delaroche 1 3_2[], ex- 
cept parameters ay, Wy and Wd are changed into 0.55 fm, 
MeV and 4.0 MeV, respectively. The proton optical poten- 
tial Up is assumed to be the same as U n - Coulomb interactions 
in the p- 209 Bi and a- 209 Bi subsystems are approximated into 
e 2 Zu Z-Qi I R, because Coulomb breakup effects are negligibly 
small in the present scattering. 
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In the pseudostate method, the discrete and discretized contin- 
uum states are obtained by diagonalizing h in a space spanned 
by £ 2 -type basis functions. As the basis function, the Gaus- 
sian or the transformed Harmonic Oscil- 
TjjGjl Eol IH is usually taken. In this 
The modelspace V is 



lator function [13 
paper, we use the Gaussian function 



then described by 



nlm I j 



(5) 




where $„/ m is the ?ith eigenstate of 6 Li with an energy e n i, a 
total spin I and its projection on the z-axis m. 

In actual calculations, the $„/ m are obtained for I n = 1 + , 
2 + and 3 + by diagonalizing h with 10 Gaussian functions for 
each coordinate in which the range parameters are taken form 
0.1 fm to 12 fm in geometric series. The $„/ m with e„j < 
20 MeV are excluded from V, since they do not affect cross 
sections of 6 Li + 209 Bi scattering. The resulting numbers of 
the discrete states are 65 (including the ground state of 6 Li), 
57 and 63 for 1 + , 2 + , and 3 + states, respectively. 

The CDCC wave function ^cdcC' w i m trie tota l an gul ar 
momentum J and its projection on the z-axis M, are ex- 
pressed as 



* JM = J2xiiL(Pm,R)/RyZ?i 



iIL 



with 



yi M = [$> n i{y c ,r c )®i L Y L {R) 



JM 



(6) 



(7) 



for the orbital angular momentum L regarding R. The 
expansion-coefficient xli where 7 = (n, I, L), describes a 
motion of 6 Li in its (n, I) state with linear and angular relative 
momenta P n i and L. Multiplying the four-body Schrodinger 
equation (0|l by y*/ M from the left and integrating it over all 
FIG. 1: Angular distribution of elastic cross section for n + 209 Bi variables except R, one can obtain a set of coupled differential 
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scattering at 5 MeV. The solid line is the result of the neutron optical 
potential t/° p . The experimental data is taken from Ref. l3lll . 



The internal Hamiltonian h of 6 Li is described by the p + 



equations for x^' 



d 2 L(L + 1) 2/x 2 
d& R~ 2 Jf U ^ R > + 



xi{Pm,R) 



n + He orthogonality condition model H33H . The Hamilto 



nian of 6 Li agrees with that of 6 He in Ref. [19], when the 
Coulomb interaction between p and 4 He is neglected. Namely, 
the Bonn-A interaction ll34ll is taken in the p-n subsystem and 
the so-called KKNN interaction f35ll is used in the p-a and 
n-a subsystems, where the KKNN interaction is determined 
from experimental data on low-energy nucleon-a scattering. 

Eigenstates of h consist of finite number of discrete states 
with negative energies and continuum states with positive en- 
ergies. In four-body CDCC, the continuum states of projectile 
are discretized into a finite number of pseudostates by either 
the pseudostate method 111 3142 lll23l425ll or the momentum-bin 
method J22ll . The Schrodinger equation ([T} is solved in a mod- 
elspace V spanned by the discrete and discretized-continuum 
states: 



2// 

ft 2 



U^(R) X ^(PnT,R) (8) 



with the coupling potentials 

Uy^R) = {yj, M \U n (R n ) + U P {R P ) + U a {R a )\yj M ), 

where \l is the reduced mass between 6 He and 209 Bi. The 
elastic and discrete breakup S'-matrix elements are obtained 
by solving Eq. (|8]l under the standard asymptotic boundary 
condition QQ. 

We also do three-body CDCC calculations by assuming a d 
+ 4 He + 209 Bi model, following Refs. JMlH. As an inter- 



(4) 



action between d and 4 He, we take the potential of Ref. (3J 
which are determined from experimental data on the ground- 
state energy ( — 1.47 MeV) and the 3 + -resonance state energy 
(0.71 MeV) of 6 Li and low-energy d-a scattering phase shifts. 
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The continuum states between d and 4 He are discretized with 
the pseudostate method B14I1 and are truncated at 20 MeV in 
the excitation energy of 6 Li from the d- 4 He threshold. The 
d- 209 Bi (type-a) optical potential (t/° p ) 1 38] is taken as U d , 
whereas U a is common between three- and four-body CDCC 
calculations. 

Results. Figure [2] shows the angular distribution of elas- 
tic cross section for 6 Li + 209 Bi scattering at 29.9 MeV. The 
dotted line shows the result of three-body CDCC calculation 
with U® p as Ud- This result, which is consistent with the pre- 
vious result of Ref. 02811 . underestimates the measured cross 
section ||26[ l27ll . The solid (dashed) line, meanwhile, stands 
for the result of four-body CDCC calculation with (with- 
out) projectile breakup effects. In CDCC calculations with- 
out 6 Li-breakup, the modelspace V is composed only of the 
6 Li ground state. The solid line reproduces the experimen- 
tal cross section, but the dashed line does not. The projectile 
breakup effects are thus significant and the present 6 Li scat- 
tering is well described by the p + n + 4 He + 209 Bi four-body 
model. This conclusion is true also for 6 Li + 209 Bi scattering 
at 32.8 MeV shown in Fig. [3] 
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FIG. 2: Angular distribution of the elastic cross section for 
6 Li + 209 Bi scattering at 29.9 MeV. The cross section is normal- 
ized by the Rutherford cross section. The dotted (dot-dashed) line 
stands for the result of three-body CDCC calculation in which ll2 F 
([/f F ) is taken as Ud- The solid (dashed) lines represent the results 
of four-body CDCC calculations with (without) breakup effects. The 
experimental data are taken from Ref. l26H27ll . 

Now we consider d-breakup in the 6 Li scattering in order to 
understand four-body dynamics of the scattering. In the limit 
of no d-breakup, the interaction between d and 209 Bi can be 
obtained by folding U n and U p with the deuteron density. This 
potential is referred to as the single-folding potential 17? F . In 
Figs. |2] and [3j the dot-dashed lines show the results of three- 
body CDCC calculations with £/f F as Ud- The results well 
simulate those of four-body CDCC calculations, i.e., the solid 
lines. This indicates that d-breakup is suppressed in the 6 Li 
scattering. Intuitive understanding of this property is as fol- 
lows. As a characteristic of the present 6 Li scattering, it is 
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FIG. 3: The same as in Fig.[2]but for 6 Li + 209 Bi scattering at 32.! 
MeV. The experimental data are taken from Ref. [26, 2^]. 
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FIG. 4: Angular distribution of the elastic cross section for d + 209 Bi 
scattering at 12.8 MeV. The solid (dashed) line stands for the result 
of three-body CDCC calculation with (without) deuteron breakup, 
whereas the dotted line is the result of the deuteron optical potential 
U2 P ■ The experimental data are taken from Ref. 13811 . 



quite peripheral in virtue of the Coulomb barrier. The scat- 
tering is dominated by the configuration in which a is located 
between d and the target, because U a is more attractive than 
Ud- In this configuration, d is out of the range of U n and 
U p , so that d-breakup is suppressed. The 6 Li elastic scatter- 
ing near the Coulomb barrier is thus well described by the 
d + a + 209 Bi three-body model, if C/| F is taken as Ud- 

Figure [4] shows the angular distribution of elastic cross sec- 
tion for d + 209 Bi scattering at 12.8 MeV. The solid and 
dashed lines stand for the results of three-body CDCC cal- 
culations with and without d-breakup, respectively, in which 
the p + n + 209 Bi model is assumed and both Coulomb and 
nuclear breakup are taken into account. The solid line re- 
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produces the data fairly well, but the dashed line does not. 
Thus d-breakup is significant for the deuteron scattering. The 
deuteron optical potential l^p p (dotted line) yields fairly good 
agreement with the data, but the radius of U2 P is larger than 
that of C/f F . This is the reason why three-body CDCC cal- 
culations with Uy F as Ud can not reproduce the measured 
elastic cross section for 6 Li + 209 Bi scattering. The difference 
between [/| F and U® p mainly comes from the fact that U® p 
includes d-breakup effects, whereas [/| F does not. 

Summary. The 6 Li + 209 Bi scattering at 29.9 MeV and 32.8 
MeV near the Coulomb barrier are well described by four- 
body CDCC based on the p + n + 4 He + 209 Bi model. This 
is the first application of four-body CDCC to 6 Li scattering. 
In the 6 Li scattering, d-breakup is strongly suppressed, sug- 



gesting that the d + 4 He + 209 Bi model becomes good, if the 
single-folding potential £/J F with no d-breakup is taken as an 
interaction between d and the target. For d+ 209 Bi scattering 
at 12.8 MeV, meanwhile, d-breakup is significant, so that the 
deuteron optical potential U2 P includes d-breakup effects. 

Four-body CDCC is applicable also for n + 6 Li scattering 
that is a key reaction in nuclear engineering. In the scattering, 
6 Li breakup into n + p + a is considered to be not negligible 
for emitted neutron spectra ll39ll . We will discuss this point in 
a forthcoming paper. 

The authors would like to thank Y. Watanabe, K. Ogata and 
K. Kat5 for helpful discussions. This work has been supported 
in part by the Grants-in-Aid for Scientific Research of Mon- 
bukagakusyou of Japan and JSPS. 



[1] S. Funada, H. Kameyama, and Y. Sakuragi, Nucl. Phys. A 575, 
93 (1994). 

[2] E. Hiyama and M. Kamimura, Nucl. Phys. A 588, c35 (1995). 
[3] T. Myo, K. Kato, S. Aoyama, and K. Ikeda, Phys. Rev. C 63, 

054313 (2001). 
[4] A. Csoto, Phys. Rev. C 49, 3035 (1994). 
[5] K. Arai, Y. Suzuki, and K. Varga, Phys. Rev. C 51, 2488 (1995). 
[6] Y. Kikuchi, N. Kurihara, A. Wano, K. Kato, T. Myo, and 

M. Takashina, Phys. Rev. C 84, 064610 (201 1). 
[7] M. Kamimura, M. Yahiro, Y. Iseri, Y. Sakuragi, H. Kameyama, 

and M. Kawai, Prog. Theor. Phys. Suppl. 89, 1 (1986). 
[8] N. Austern, Y Iseri, M. Kamimura, M. Kawai, G. Rawitscher, 

and M. Yahiro, Phys. Rep. 154, 125 (1987). 
[9] M. Yahiro, K. Ogata, T. Matsumoto, and K. Minomo, 

larXiv: 1203.53921 [nucl-th] . 
[10] N. Austern, M. Yahiro, and M. Kawai, Phys. Rev. Lett. 64, 2649 

(1989). 

[11] N. Austern, M. Kawai, and M. Yahiro, Phys. Rev. C 53, 314 
(1996). 

[12] A. Deltuva, A. Moro, E. Cravo, F. Nunes, and A. Fonseca, Phys. 

Rev. C 76, 064602 (2007). 
[13] A. M. Moro, J. M. Arias, J. Gomez-Camacho, I. Martel, 

F. Perez-Bernal, R. Crespo, and F. Nunes, Phys. Rev. C 65, 

011602(R) (2001). 
[14] T. Matsumoto, T. Kamizato, K. Ogata, Y Iseri, E. Hiyama, 

M. Kamimura, and M. Yahiro, Phys. Rev. C 68, 064607 (2003). 
[15] T. Matsumoto, E. Hiyama, K. Ogata, Y. Iseri, M. Kamimura, 

S. Chiba, and M. Yahiro, Phys. Rev. C 70, 061601(R) (2004). 
[16] T. Egami, K. Ogata, T. Matsumoto, Y. Iseri, M. Kamimura, and 

M. Yahiro, Phys. Rev. C 70, 047604 (2004). 
[17] M. Rodriguez-Gallardo, J. M. Arias, J. Gomez-Camacho, A. M. 

Moro, I, J. Thompson, and J. A. Tostevin, Phys. Rev. C 72, 

024007 (2005). 

[18] A. M. Moro, F. Perez-Bernal, J. M. Arias, and J. Gomez- 
Camacho, Phys. Rev. C 73, 044612 (2006). 

[19] T. Matsumoto, T. Egami, K. Ogata, Y. Iseri, M. Kamimura, and 
M. Yahiro, Phys. Rev. C 73, 051602(R) (2006). 

[20] M. Rodriguez-Gallardo, J. M. Arias, J. Gomez-Camacho, R. C. 



Johnson, A. M. Moro, I. J. Thompson, and J. A. Tostevin, Phys. 

Rev. C 77, 064609 (2008). 
[21] A. M. Moro, J. M. Arias, J. Gomez-Camacho, and F. Perez- 
Bernal, Phys. Rev. C 80, 054605 (2009). 
[22] M. Rodriguez-Gallardo, J. M. Arias, J. Gomez-Camacho, A. M. 

Moro, I. J. Thompson, and J. A. Tostevin, Phys. Rev. C 80, 

051601(R) (2009). 
[23] T. Egami, T. Matsumoto, K. Ogata, and M. Yahiro, Prog. Theor. 

Phys. 121, 789 (2009). 
[24] T. Matsumoto, T. Egami, K. Ogata, and M. Yahiro, Prog. Theor. 

Phys. 121, 885 (2009). 
[25] T. Matsumoto, K. Kato, and M. Yahiro, Phys. Rev. C 82, 

051602(R) (2010). 
[26] E. F. Aguilera etal, Phys. Rev. Lett. 84, 5058 (2000). 
[27] E. F. Aguilera etal, Phys. Rev. C 63, 061603 (2001). 
[28] N. Keeley, J. M. Cook, K. W. Kemper, B. T. Roeder, W. D. 

Weintraub, F. Marechal, and K. Rusek, Phys. Rev. C 68, 054601 

(2003). 

[29] K. Rusek, I. Martel, J. Gomez-Camacho, A. M. Moro, and 

R. Raabe, Phys. Rev. C 72, 037603 (2005). 
[30] A. R. Barnett and J. S. Lilley, Phys. Rev. C 9, 2010 (1974). 
[31] J. Annand, R. Finlay, and P. Dietrich, Nuclear Physics A 443, 

249 (1985). 

[32] A. J. Koning and J. P. Delaroche, Nucl. Phys. A 713, 231 
(2003). 

[33] S. Saito, Prog. THeor. Phys. 41, 705 (1969). 

[34] R. Machleidt, Adv. Nucl. Phys. 19, 189 (1989). 

[35] H. Kanada, T. Kaneko, S. Nagata, and M. Nomoto, Prog. Theor. 

Phys. 61, 1327 (1979). 
[36] R. A. D. Piyadasa, M. Yahiro, M. Kamimura, and M. Kawai, 

Progress of Theoretical Physics 81, 910 (1989). 
[37] Y. Sakuragi, M. Yahiro, and M. Kamimura. 
[38] A. Budzanowski, L. Freindl, K. Grotowski, M. Rzeszutko, 

M. Slapa, J. Szmider, and P. Hodgson, Nuclear Physics 49, 

144 (1963). 

[39] T. Matsumoto, D. Ichinkhorloo, Y. Hirabayashi, K. Kato, and 
S. Chiba, Phys. Rev. C 83, 064611 (2011). 



